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ABSTRACT 

We compare two different methods of constraining the characteristic velocity and spatial 
scales of gas motions in the X-ray bright, nearby Centaurus cluster, using new deep (760ks) 
Chandra observations. The power spectrum of excess surface brightness fluctuations in the 
0.5-6.0 keV band in a sector to the west is measured and compared to theoretical expectations 
for Kolmogorov index fluctuations. The observed power spectrum is flatter than these expec¬ 
tations, and the surface brightness fluctuations are around the 8 percent level on length scales 
of 2 kpc. We convert the 2D power spectrum of fluctuations into a 3D power spectrum using 
the method of Churazov et ah, and then convert this into constraints on the one-component 
velocity of the gas motions as a function of their length scale. We find one-component veloc¬ 
ities in the range 100-150 km/s on spatial scales of 4-10 kpc. An independent constraint on 
the characteristic velocity and length scales of the gas motions is then found by considering 
the diffusion coefficient needed to explain the distribution of metals in the Centaurus cluster, 
combined with the need to balance the rate of gas cooling with the rate of heat dissipated by 
the gas motions. We find that these two methods of constraining the velocity and length scales 
of the gas motions are in good agreement. 

Key words: galaxies: clusters: individual: Centaurus - galaxies: clusters: intracluster medium 
- intergalactic medium - X-rays: galaxies: clusters. 


1 INTRODUCTION 

In the centres of galaxy clusters the cooling time of gas is much 
shorter than the Hubble time. Feedback from the central AGN is 
needed to counteract this cooling to prevent the formation of a cool¬ 
ing flow and limit the rate of star formation. The precise way in 
which the AGN feedback energy is dissipated into the intracluster 
medium (ICM), and how this can so accurately balance the rate of 
cooling, remains under debate. 

Turbulence is one mechanism through which the energy from 
the central AGN could be dissipated into the ICM to counteract 
cooling. Direct constraints on the level of turbulence in the cores 
of galaxy clusters remain challenging with existing X-ray instru¬ 
ments. Measuring the width of emission lines seen with the RGS 
(reflection grating spectrometer) instrument on XMM- Newton can 
provid e upper limits on the level of gas motions. I Sanders et al.l 
yOlOj) used this method to find a limit of 274 km s _1 for ran¬ 
dom gas motions in the core of Abell 1835, while in the sample of 
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clusters studied in I Sanders et al.l l l201 lh 15 clusters were found for 
which the upper limits gave a turbulence to thermal energy density 
ratio of less than 20 percent. 

Recent high spatial resolution studies of nearby clusters have 
investigated spatial fluctuations in the X-ray emission from the cen¬ 
tres of galaxy clusters, which can place limits on the level of gas 
motions in the ICM due to turbulence. IChurazov et alj < t2012h stud¬ 
ied the X-ray surface brightness fluctuations in the core of the non¬ 
cool core Coma cluster with C handra and XMM-Newton data, us¬ 
ing the A-variance technique istutzki et al.ll 1 998i) to characterise 
the magnitude of fluctuations as a function of their spatial scale. 
They found that the characteristic amplitude of excess density fluc¬ 
tuations varied from ~5 percent on scales of ~30 kpc, up to 7-10 
per cent on much larger scales of ~ 500 kpc. 

ISanders & Fabian! J2012I ) investigated the X-ray surface 
brightness fluctuations in the cool-core galaxy cluster AWM7, and 
found implied 3D density density fluctuations with a standard de¬ 
viation of around 4 percent, with 3D pressure variations of at most 
4 per cent as well. They compared the shape of the power spectrum 
of excess surface brightness fluctuations with that expected from 
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Kolmogorov index fluctuations, and found the observed spectrum 
to be flatter than these expectations. By comparing the northern and 
southern sectors of AWM7, they found significant spatial variations 
in the fluctuations. In the northern sector studied, there was a factor 

of 2 less A-vari ance relative to th e southern secto r studied. _ 

_ Recently, IZhuravleva et al.l d2014l> and Izhuravleva et al.l 

(120151) have studied the X-ray s urface brightness fluctua tions in 
the Perseus and Virgo clusters. In Izhuravleva et ail d2014l) the 2D 
power spectrum of surface brightness fluctuations is converted into 
a 3D power spectrum, which is subsequently converted into the 
one-component velocities of the gas motions needed to cause the 
fluctuations, which are found as a function of their spatial size. This 
method was applied to annuli of increasing radius, and the charac¬ 
teristic one-component velocity was found to decrease from around 
200 km/s in the central 30 kpc to 100 km/s at a distance of 200 kpc 
from the centre. They then calculated the turbulent heating rate ex¬ 
pected for gas motions with the velocity and spatial sizes found. For 
the Perseus and Virgo clusters, they find that the rate of turbulent 
heating deduced using this method is roughly equal to the cooling 
rate of the gas in all of the radial annuli studied, adding support 
to the idea that turbulence plays a role in balancing gas cooling in 
cluster cores. 

Independent constraints on the velocity and spatial scales of 
the gas motions can also be found from the distribution of metals 
in the abundance peak around the brighte st cluster galaxy, BCG 
dRebusco et al.ll2005l . iRebusco et alj|2006h . This method assumes 
that, as the metals are continuously deposited into the ICM from 
stars in the BCG over the life time of the cluster, the profile of 
the metal abundance from this deposition follows the light profile 
of the BCG. The metals then diffuse outwards to form the much 
broader metal abundance peaks in cluster cores that we see today. 
By assuming that metal diffusion has proceeded through stochas¬ 
tic gas motions, this allows the metal diffusion coefficient, D to 
be constrained, which is proportional to the product of the charac¬ 
teristic velocity, v, and length scales, l, of these gas motions. If it 
is then assumed that the heating rate from the dissipation of these 
stochastic gas motions balances the gas cooling rate, v and l can be 
constrained. IRebusco et al.l2005l used this technique for the Perseus 
cluster and found that a diffusion coefficient of 2 x 10 29 cm 2 s -1 
is needed to explain the observed abundance profile. They deduced 
that the turbulent gas motions in Perseus have velocities of ~300 
km s' 1 on spatial scales of 20 kpc. 

Here we investigate new, deep (760 ks) Chandra observations 
of the cool-core Centaurus cluster (PI J.S. Sanders) to compare the 
constraints on the characteristic velocity and spatial scales of gas 
motions obtained by considering the X-ray surface brightness fluc¬ 
tuations and the distribution of metals in the ICM. A detailed study 
of these new Chandra observations will be presented in Sanders 
et al. (submitted). Centaurus’ proximity (it is the second nearest 
galaxy cluster), combined with its high X-ray flux makes it an ideal 
target for exploring X-ray surface brightness fluctuations. 

iGraham et al.l (2006b have previously used the existing 200ks 
Chandra observ ation tojnvestigate metal diffusion following the 
method of IRebusco et al.l 12003 ). The azimuthally averaged metal 
abundance profile for Centaurus has a prominent drop between 20- 
50 kpc from the core which requires a comparatively low diffusion 
coefficient to maintain. IGraham et al.1 d2006b found a best fitting 
diffusion coefficient in the range 2-4 xlO 28 cm 2 s -1 , nearly an or¬ 
der of magnitude lower than that found for Perseus in lRebusco et al .1 
(2005). The very high meta llicity around the BCG of Centaurus, 
which peaks at nearly 3 Zq jFabian et al.l2005l ). makes it an excel¬ 
lent target for applying this method. 



Figure 1. Exposure corrected 0.5-6.0 keV mosaic image of the Centaurus 
cluster. The white sector region shown is used for the analysis of the sur¬ 
face brightness fluctuations. Point sources have not been removed from this 
image. 


Table 1. Chandra observations used 


Obsid 

Cleaned exposure/ ks 

Date 

504 

26 

2000-05-22 

505 

10 

2000-06-08 

4190 

34.1 

2003-04-18 

4191 

33.8 

2003-04-18 

4954 

85.4 

2004-04-01 

4955 

43.1 

2004-04-02 

5310 

48.7 

2004-04-04 

16223 

176.2 

2014-05-26 

16224 

41.2 

2014-04-09 

16225 

29.7 

2014-04-26 

16534 

55 

2014-06-05 

16607 

44.6 

2014-04-12 

16608 

33.3 

2014-04-07 

16609 

81.7 

2014-05-04 

16610 

17.1 

2014-04-27 

Total 

760 



We use a standard ACDM cosmology with Ho = 70 km s -1 
Mpc” 1 , flju = 0.3, Da=0.7. For the redshift of the Centaurus 
cluster, 2 = 0.011, the luminosity distance is 49.5 Mpc and the 
angular scale is 0.235 kpc/arcsec. All errors unless otherwise stated 
are at the 1 a level. 


2 DATA 

The new deep Chandra observations of Centaurus will be presented 
and explored in depth in Sanders et al. (submitted). The observa¬ 
tions used are tabulated in table[T] The data were reprocessed using 
CIAO v4.6. Lightcurves for each observation were examined, with 
any periods deviating from the standard deviation by more than 2.5 
a removed, leaving a total clean exposure of 760ks. 
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Figure 2. Left: Actual 0.5-6.0keV band image from the sector region shown in Fig. [T] with point sources removed. Right: Simulated 0.5-6.0keV band image 
from the same sector, obtained by modelling the surface brightness of this region with 10 ellipsoids. 


3 ANALYSIS 

The exposure corrected and background subtracted 0.5-6.0 keV 
band stacked image of Centaurus is shown in Fig. [T] The stacked, 
exposure corrected image was produced following the methods de¬ 
scribed in lWalker et al .1 a201 4 ) . To study the surface brightness fluc¬ 
tuations we investigate the region shown by the white sector, which 
spans a radial range of 23-43 kpc from the centre of the BCG. This 
region lies between the core and the large western cold front. It has 
a high X-ray surface brightness over a large area and is relatively 
smooth in its X-ray emission. It was chosen to avoid any surface 
brightness features which are obviously due to feedback from the 
core (i.e. the cavities and filaments) or due to gas sloshing (i.e. the 
eastern and western cold fronts and the ‘bay’ feature to the south 
west of the core). 

The left hand panel of Fig.[2]shows the sector region by itself, 
with point sources removed. Point sources were identified by the 
tool wavdetect, using a range of wavelet sizes between 1 and 16 
pixels. The large scale surface brightness distribution in this region 
was modelled by fitting ten ellipsoids to it, producing a smooth 
model image. This ten ellipsoid model was used because it was 
found to be the simplest way of accurately modelling the large scale 
surface brightness distribution of the region. The power spectrum 
of this model on large scales matches that of the real observations 
very closely. This smooth model image was multiplied by the ex¬ 
posure map, and then a simulated image in the 0.5-6.0 keV band 
was obtained by making a Poisson realisation. This was then ex¬ 
posure corrected by dividing by the exposure map, producing the 
simulated image shown in the right hand panel of Fig. [2] 


In Fig. [3] we show the ratio of the real image (left) and the 
simulated image (right) to the smooth image we obtained by fitting 
10 ellipsoids to this region. There are clear structures in the residu¬ 
als of the real data to this smooth model (left hand panel) at the 10 
percent level. 


4 COMPARING FLUCTUATION POWER SPECTRA 

To compare the shape and magnitude of the observed power spec¬ 
trum of surface brightness fluctuations with theoretical expecta¬ 
tions for hydrodynamical turbulence, we follow the approach of 
ISanders & Fabianl J2012h . We use the simulated image of the sec¬ 
tor, obtained by fitting 10 ellipsoids to the region, and add to it 
surface brightness fluctuations with a powerlaw power-spectrum. 
These surface brightness fluctuations are generated by finding the 
fast Fourier transform of a cube of dimensions 700x700x700, 
which is larger than the size of the region of the Centaurus image 
explored (in terms of pixels). Each pixel in the cube was scaled by 
a powerlaw power spectrum of index -11/3 (the shape expected for 
the electron density power spectrum for a Kolmogorov process), 
using a normal distribution for the real and imaginary components. 
Following the inverse transform, the real component for each pixel 
was taken. The powerspectrum was truncated to remove structure 
on scales larger than 30 kpc (the radial range spanned by the sec¬ 
tor examined). We ensured that this truncation scale is consistent 
with the use of the 10 ellipsoids model for removing large-scale 
perturbations. 

We deproject the surface brightness profile from the Cen- 
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Figure 3. Left: Ratio of the actual data to the best fitting 10 ellipsoidal model. Right: Ratio of the simulated data to the best fitting 10 ellipsoidal model. 


taurus observations to obtain the deprojected emissivity model, 
E(x, y, z). We then integrate the cubic fluctuation, F(x, y, z) (nor¬ 
malised to have a standard deviation of 1), field along the 2 -axis, 
multiplying by the deprojected emissivity profile and a scale factor 
determining the magnitude of input fluctuations, e, 

S(x,y) = J[l + eF(x,y,z)]E(x,y,z)dz (1) 

from which we obtain a simulated image, S(x,y) containing 
powerlaw fluctuations of magnitude 4, 8 and 16 percent in surface 
brightness. The Chandra exposure maps were then used to create 
Poisson realisations of these images, allowing us to compare these 
to the actual observations. 

We then used the delta variance technique to find the power 
spectra of the simulated images to which fluctuations of magnitude 
4, 8 and 16 percent in surface brightness have been added. These 
magnitudes correspond to the RMS of the surface brightness over 
the image. These power spectra are shown in Fig. [4] The power 
spectrum of the simulated ellipsoidal model with zero fluctuations 
has been subtracted from these power spectra, allowing us to see 
the excess power that is just due to the fluctuations we have added. 
The error ranges of the power spectra were calculated by simulat¬ 
ing 1000 Poisson realisations of each image and finding the stan¬ 
dard deviation of the power spectra from these. We performed tests 
to ensure the results are not sensitive to variations in the procedure 


used. The analysis was repeated but this time the model was sub¬ 
tracted from the data first, and then the power spectrum was found. 
The results were found to be the same as those obtained by sub¬ 
tracting the power spectrum of the model from that of the data. 

The power spectrum of the actual data (shown as the cyan re¬ 
gion in Fig. 0 is then compared to these simulations of different 
fluctuation levels. We see that the power spectrum of the real data 
is flatter in shape than the spectra made assuming fluctuations using 
a standard Kolmogo rov powerlaw index of -1 1/3. This is similar to 
what was found by ISanders & Fabiar] < 120121 ) for AWM7. It may 
be possible that unresolved compact sources contribute to the ap¬ 
parent flatness of the observed spectrum. The 10 ellipsoidal model 
used only has a noticeable affect on scales larger than around 10 
kpc, so it is unlikely that the use of this model contributes to the 
observed flattened shape. 


5 VELOCITY POWER SPECTRA 

Following the method detailed in IChurazov et al .1 d2012h . we can 
convert the 2D power spectrum obtained from our image analysis, 
P 20 , into a 3D power spectrum P 3 d using 

P2o(k) «4P 3D (fc) J\W{k z )\ 2 dk z , (2) 


© 0000 RAS, MNRAS 000, 000-000 















Centaurus gas motions 5 


10"' 8 



I 


i cr 20 


0.5 1 2 3456789 

Length/kpc 



Figure 4. The power spectrum for the excess surface brightness fluctuations in the real data, from which the spectrum of the smooth model has been subtracted, 
is shown as the cyan region. This is compared to the power spectra of simulated images to which 4, 8 and 16 percent Kolmogorov index fluctuations have been 
added. 
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Figure 5. Left: Constraints on the one component velocity of the gas motions as a function of their wavenumber, derived from the power spectrum of excess 
surface brightness fluctuations. The solid cyan region shows the region where systematic and statistic uncertainties least affect the measurement. Right: The 
cyan region shows the turbulent heating rate inferred from the surface brightness fluctuations against the cooling rate of the sector region studied. This agrees 
reasonably well with equilibrium between heating and cooling, which is shown by the dashed black line. 


where | W(k z ) | 2 is the ID power spectrum of the normalised emis- 
sivity, rj(z), along the line of sight 

n(x,y,z) = 

and for the small area of the cluster we consider the dependence on 

x and y can be neglected, such that r/(x, y , z) sa r/(z). _ 

We can then follow the approach of IZhuravleva et ahl d2014l) 
and lZhuravleva et all d20 1 5l ) to convert this 3D power spectrum into 
where 7 /( 2 ) is related to the density distribution n^x, y, z) by the one-component velocity, Vi } k power spectrum. For the strati- 


W(k z ) — J rj(z)e z2nzkz dz. 


(3) 


nl{x,y,z) _ nl{x,y,z) 
f n o(x, V, z')dz' Io(x,y) 


(4) 
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fied atmospheres of relaxed galaxy clusters, and within the inertial 
range of scales, the one-component velocity is expected to be pro¬ 
portional to the amplitude of density fluctuations, Azo = at 
each length scale l = 1/k, so that we can write 

A 3D = 5 -^=ni 1 -^ ( 5 ) 

p 0 Cs 

where c s is the sound speed. [Zhuravleva et ali d2014t) found using 
cosmological simulations of galaxy clusters that the average value 
of the proportionality constant r/ 1 is consistent with unity (r/i = 
1 ± 0.3). 

The amplitude of the density fluctuations is related to the 3D 
power spectrum by 


A 3 d = \J PzD(k)4nk 3 (6) 

which allows us to find the amplitude of the one-component ve¬ 
locity of the gas motions as a function of wavenumber (fc = 1/1) 
from 


Vi ,fc « c s \/P zd (fc) 47rfc 3 /r)i (7) 

The resulting plot of Vi t k against wavenumber is shown in the 
left hand panel of Fig. [5] where the solid cyan region shows where 
the constraint is least affected by systematic and statistic uncertain¬ 
ties. 

In the right hand panel of Fig. [5] we show the turbulent heat¬ 
ing rate for the sector derived from the surface brightness fluctu¬ 
ations using equation [TO] plotted against the cooling rate for the 
same region. The two appear to agree reasonably well. This is sim- 
ilar to the findings for the Perseus and Virgo clusters studied in 
Zhuravleva et al.U2014h . 


iGraham et alj 120061) found that the best fitting diffusion coef¬ 
ficient had a radial dependence, decreasing with increasing radius. 
They found that D could be described by 


D 


4 x 10 28 


n H (r) 

n H (ro) 


2 -1 
cm s 


(9) 


where uh is the hydrogen density profile, and r 0 = 25kpc. Over 
the radial range of the sector we consider (23-43 kpc), D decreases 
from 4.7 to 1.5xl0 28 cm 2 s -1 . We can then plot the constraints on 
this relation between v and l , and this is shown with the black lines 
in the left hand panel of Fig. [6] 

Another relation between v and l can be obtained by requiring 
the heat diss ipation rate^of the gas motions to balance the cooling 
rate. From lRebusco et alj fl2005l) . the heating rate can be written as 


r diS s ~ OApv 3 /i (10) 

from which we obtain 

n(r) 2 A(T(r), A(r)) ~ OApv 3 /l (11) 

The red lines in the left hand panel of Fig. [6] show the constraints 
on v and l from this relation, where the range takes into account 
the range of gas density, temperature and metal abundance over the 
radial range of the sector considered. The constraints from equa¬ 
tions [8] and QT] overlap in a region where they are both consistent, 
with heating balancing cooling and the metal diffusion coefficient 
agreeing with observations. 

In the right hand panel of Fig. [6] we compare these constraints 
on v and l with those obtained earlier in Fig. [5] by considering the 
power spectrum of excess surface brightness fluctuations. The two 
different methods appear to be in good agreement with one another. 


6 CONSTRAINTS ON VELOCITY AND SPATIAL 
SCALES CONSIDERING METAL DIFFUSION 
THERMAL EQUILIBRIUM 

Independent constraints on the velocities and spatial scales of the 
gas motions can be obtained by considering the metal diffusion co¬ 
efficient, D, and the level of gas m otions needed to offset cool¬ 
ing in the ICM. Following iRebusco et alj J20051) (who considered 
the Perseus cluster), under the assumption that metal diffusion is 
caused by stochastic gas motions on a length scale smaller than 
that distance from the centre of the cluster, r, we can write the dif¬ 
fusion coefficient in terms of the length scale, l and velocity scale, 
v of the gas motions 


D ~ O.lluZ 


(8) 


IGraham et alj d2006lf applied the stoc hastic model of metal 
transport developed by Rebusco et al J <120057 to the Centaurus clus¬ 
ter to obtain constraints on D. This method constrains D by using 
the central abundance peak in galaxy clusters. The central abun¬ 
dance peak is likely produced by the stars in the BCG following the 
formation of the cluster, so the initial distribution of metals (when 
deposited) would be expected to follow the distribution of stars in 
the BCG. The observed metal peak is much broader than the light 
profile of the BCG, and if we assume the metals have been trans¬ 
ported outwards by a stochastic process, constraints on D can be 
made by comparing the two profiles. 


7 CONCLUSIONS 

We have examined the A-variance power spectrum of the excess 
X-ray surface brightness fluctuations in a large sector to the west of 
the core of the Centaurus cluster, spanning a radial range between 
23-43 kpc from the center of the cluster. When compared to the 
expected shape of the spectrum from Kolmogorov index fluctua¬ 
tions, the obs erved spectrum app e ars fla tter. This is similar to what 
was found in iSanders & Fabianl d2012l) for the fluctuation power 
spectrum of AWM7. The magnitude of the surface brightness fluc¬ 
tuations is roughly 8 percent on scales of 2 kpc. 

We have converted the 2D power spectrum into a 3D power 
spectrum and then found constraints on the characteristic velocity 
of the gas mot ions as a function of t heir spatial scaling fo l lowin g 
the methods of lChurazov et al.l J2012h and lZhuravleva et ahl J 201 4t) . 
On scales between 4-10 kpc where the method is least sensitive to 
systematic and statistical errors, the one-component velocity of the 
gas motions is calculated to be in the range 100-150 km/s using 
this method. The turbulent heating rate from such gas motions is 
consistent with what is needed to balance the cooling rate, as wa s 
found for the Perseus and Virgo clusters in lZhuravleva et alj 1 120141) . 

Independent constraints on the characteristic velocity and size 
scale of the gas motions have been obtained by measuring the dif¬ 
fusion coefficient needed to explain the present day distribution of 
metals, and invoking equilibrium between the turbulent dissipation 
rate and the gas cooling rate. This involves assuming that when the 
metals were first deposited by the stars in the BCG their abundance 
profile followed the light profile of the BCG. The constraints on v 
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Figure 6. Left: Th e sol i d blac k line shows the relation between v and l obtained by fixing the metal diffusion coefficient in equation [8]to the best fitting value 
found in iGraham et alj J20061) for a radius at the midpoint of the sector. The dashed black lines show the allowed range for this relation when we take into 
account the the radial variation of the D in the sector studied, and an uncertainty of a factor of 2 in the constant coefficient in equation [8] The solid red line 
shows the relation obtained by equating the turbulent heating rate with the cooling rate in equation II11 The dashed red lines show the range of this relation 
when the radial dependence of the gas density and cooling function are factored in, along with an uncertainty of a factor of 3 in the constant coefficient in 
equation EL The region where these two methods overlap provides a constraint on v and l for the gas motions Right: Here we overplot the constraints on v 
and l from the left panel onto those obtained obtained in Fig.[5]by considering the power spectrum of excess surface brightness fluctuations. 


and l obtained using this method are in reasonable agreement with 
those obtained from the X-ray surface brightness fluctuations. 

It is possible that some fraction of the observed surface bright¬ 
ness fluctuations arise from sound waves rather than turbulence 
(Sanders et al. submitted). We investigate this in Appendix A, in 
which we compare the power spectrum from the whole sector with 
the power spectrum obtained by removing any features resembling 
sound waves. 
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APPENDIX A: LINEAR FEATURES 

Sanders et al. (submitted) have explored applying a Gaussian gradi¬ 
ent magnitude (GGM) filter (similar to a Sobel filter) to the Chandra 
X-ray images, which emphasises gradients in the surface bright¬ 
ness. Sanders et al. (submitted) have identified several linear fea¬ 
tures in the northern part of the sector we have explored, which may 
be sound waves. This method is detailed in depth in in Sanders et 
al. (submitted). To investigate the affect of these, we repeated our 
analysis by excluding these features by examining just the southern 
half of the sector. The results are shown in Fig. lAll In the left hand 
panel we compare the power spectrum of the whole sector (cyan 
region) with that of the southern half of the sector, which excludes 
the linear features (green region). We see that there is slightly less 
power in the southern half of the sector, though given the errors 
its power spectrum is consistent with that of the whole sector. This 
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Figure Al. Le/bComparing the power spectrum of the surface brightness fluctuations in the southern half of the sector away from any linear features (green 
line shaded region) with the power spectrum of the whole sector (cyan region). Right:C omparing the constraints on the one component velocity of the gas 
motions as a function of their length scale for the two regions. The black and red dashed lines show the constraints shown earlier in Fig. [6] obtained by 
considering metal diffusion and the heating-cooling balance. 


reduction in power lowers the velocities of the gas motions calcu¬ 
lated from the surface brightness fluctuations, as shown in the right 
hand panel comparing the velocity constraints for the southern half 
of the sector with the whole sector. Given the errors the two are 
still consistent with one another, and also still consistent with the 
constraints on the velocity and length scales obtained by consider¬ 
ing metal diffusion and the balance between heating and cooling 
(dashed black and red lines respectively). 
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